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ABSTRACT 

The present work purports to identify candidate carriers of the UIBs. This requires 
a procedure for the computation of the emission spectrum of any given candidate. The 
procedure used here consists in exciting the carrier into a state of internal vibration, 
waiting until the system has reached dynamic equilibrium and, then, monitoring the 
time variations of the overall electric dipole moment associated with this vibration. 
The emission spectrum is shown to be simply related to the FT of these variations. 
This procedure was applied to more than 100 different chemical structures, inspired 
by the exhaustive experimental and theoretical analyses of Kerogens, the terrestrial 
sedimentary matter, which is known to be mainly composed of C, H, O, N and S 
atoms. From this data base, 21 structures were extracted, which fall in 4 classes, each 
of which contributes preferentially to one of the main UIBs. Summing their adequately 
weighted spectra delivers an emission spectrum which indeed exhibits the main UIB 
features (allowing for computational errors inherent in the chemical simulation code). 
By changing the weights, it is possible to change the relative band intensities so as to 
mimic the corresponding (moderate) changes observed in the sky. The defects of the 
present simulation are discussed, and directions for improvement are explored. 
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1 INTRODUCTION 

The quest for a realistic model of the structure, composition 
and excitation of the carriers of UIBs (Unidentified Infrared 
Bands) has been going on for more than three decades now. 
An essential requirement is, for such a model, to correctly 
reproduce the spectral features of UIB emission, which, by 
now, are very well documented in a large variety of astro- 
nomical sites. Among the latest publications, note the ex- 
haustive analysis of a rich collection of spectra from star- 
forming regions (Smith et al. 2007). Despite some variations 
from object to object, there are enough similarities that it 
has been possible to partition most of the observed spectra 
into a small number of classes (see, for instance, Peeters et 
al. 2002, 2004a; van Diedenhoven et al. 2004), thus allow- 
ing systematic comparisons to be made with laboratory or 
theoretical spectra. 

Most present day candidate model carriers are big 
molecules of widely different sizes (see Kwok and Sand- 
ford 2008, Tielens 2008, Kwok 2009), the size being limited, 
of course, by experimental or computational constraints, or 
even by theoretical considerations. While the experimental 
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production and spectral analysis of single macromolecules 
have proven to be daunting, modern computing capabilities 
make it easier to use theory or chemical simulation for mod- 
eling purposes. 

The most usual computational way of determining the 
IR (InfraRed) spectrum of a given structure (molecule or 
grain) is to perform a Normal Mode Analysis (NMA; see 
Wilson et al. 1955). This delivers all the different modes of 
atomic vibration of infinitesimal amplitude, characteristic of 
the structure. For a grain made up of N chemically bound 
atoms, there are 3N-6 such modes, each characterized by 
its frequency and the directions and velocities of the excur- 
sion of each atom from its equilibrium position along the 3 
space coordinates. Further treatment yields the mode ir in- 
tensity (or line intensity, or integrated band intensity), which 
is proportional to the absorbance at the corresponding fre- 
quency, the quantity commonly measured in the laboratory. 
NMA can be performed using chemical simulation codes 
based on different theories and approximations: Molecular 
Mechanics (MM+, Ambers, BIO+, etc.) and quantum chem- 
ical codes: Semi-empirical (AMI, PM3, etc.), ab initio (e.g. 
DFT) methods, which differ widely in accuracy and time 
consumption. These are available with commercial packages 
such as HyperChem, Spartan, Gaussian, etc. 



2 R. Papoular 



As early as 1992, Brenner and Barker (|1992|) warned 
that "unambiguous identification of the IS (InterStellar) 
emitter cannot be established just on the basis of matching 
the emission frequencies to laboratory absorption spectra". 
For emission spectra differ in many respects from absorp- 
tion spectra: the features of the former are generally shifted 
and broadened relative to those of the latter, because of 
anhormonicity effects, which set in as soon as the grain is 
excited; cor relatively, they depend on the excitation process 
(photon absorption, chemical reaction, etc.). Brenner and 
Barker themselves computed the emission of excited ben- 
zene and naphtalene based on a) the thermal equilibrium as- 
sumption (statistical distribution of energy among molecular 
states), b) Einstein's spontaneous emission coefficients, as 
deduced from measured (relative) absorption coefficients at 
a few characteristic spectral frequencies of the two molecules 
(see also Allamandola et al. 1989). Intrinsic bandwidth and 
underlying emission continuum were assumed. Several off- 
springs of this procedure have been used (e.g. Verstraete et 
al. 2001, Boersma et al. 2010, 2011); the state of the art 
in this domain is clearly summarized by Bauschlicher et al. 
(|20TU)) . 

Addressing the same issue, Papoular (2001 J, (2012) pro- 
posed a distinct procedure which bypasses NMA altogether. 
Based on the unversal relation between emission, electric 
dipole moment and frequency, he showed that the emission 
spectrum of an excited molecule is simply related to the 
energy content and overall dipole moment of the structure. 
The proposed procedure hence consists in monitoring the 
variations in time of the energy content and overall dipole 
moment of a selected model structure after it has been ex- 
cited. The FT (Fourier Transform) of these variations de- 
liver the emission spectrum. As the specific anharmonicity 
of the molecule, and all its effects, are reflected in its in- 
ternal motions, it also affects both the energy and dipo- 
lar spectra, and, hence, the emission spectrum. As a conse- 
quence, the absorption line spectrum of NMA is replaced by 
an emission band spectrum. No assumptions are added to 
shift or broaden the bands artificially. Plateaus arise natu- 
rally from overlapping band wings. If the excitation energy 
is high enough, overtone and combination bands automati- 
cally emerge from the analysis. 

The anharmonicity is generally built in the dynamics 
code used for chemical simulation. In semi-empirical codes, 
it is generated by a system of parameters tailored by analogy 
with measured properties of structures made of the partic- 
ular atomic species and bonds of interest. The PM3 code 
used here is best fit to simulate hydrocarbon molecules. 
Semi-empirical and Molecular Mechanics commercial codes 
were used to demonstrate relaxation after excitation (see 
Papoular 2001 and 2002), as well as anharmonicity effects: 
line shift, mode- locking, combination frequencies, Fermi- 
resonance (see Papoular 2006) and line broadening (Pa- 
poular 2012). 

Both DFT (Density Functional Theory) and Semi- 
empirical quantum chemistry methods deliver accurate Nor- 
mal Mode intensities and frequencies. However, DFT meth- 
ods are basically designed to deal with equilibrium states, 
so they are not directly applicable to the analysis of dipole 
moment variations. To my knowledge, even Time-Dependent 
DFT (TDDFT) has been applied only to adiabatic transi- 
tions between equilibrium states, not to dynamics. 



Although, the procedure described above provides a 
sound theoretical transition from NMA to emission spec- 
tra and yields a degree of line broadening, this does not 
allow to mimic the wide mid-IR bands observed in the sky, 
unless very large aggregates of atoms are considered, or the 
emission spectra of a very large number of different medium- 
sized molecules are summed up. The former option is exces- 
sively computer-time consuming, and is not pursued fur- 
ther here. As for the latter, it is preferably started with a 
survey of the Normal Mode spectra (easier to compute!) of 
structures built up with the most common chemically ac- 
tive atoms: C, H, O, N and S (CHONS for short). Even 
with this restriction, the number of possible combinations is 
daunting. Papoular (|2010p therefore sought cues in kerogens 
and coals, whose spectra have features in common with the 
UIBs, and which have been abundantly modeled using labo- 
ratory analysis and chemical codes (e.g. Durand 1980, Behar 
and Vandenbroucke 1986, Carlson 1992, Speight 1994). It 
turned out that some types of structures have spectral lines 
falling mostly within one or two UIBs. Changing slightly 
the structure or composition slightly displaces the spectral 
lines, thus helping to build a model band. More than a hun- 
dred normal mode spectra were computed and distributed 
over eight families of structures, each contributing mainly to 
one or two spectral regions. By combining these spectra in 
the right proportions, it was possible to synthesize spectra 
bearing some resemblance with UIB spectra. 

These computational experiments led to the selection 
of a number of structural types exhibiting spectral features 
which appeared to be promising for our present purposes. 
In the present work, the emission spectra of many more 
structures of these types were computed according to the 
procedure outlined above. Each was given a weight and their 
weighted sum compared with UIB spectra. The details and 
outcome of this endeavour are given below. 

Section 2 lays down the theoretical basis of the proposed 
method of spectral emission calculation; it is complemented 
by Appendix A. The details of its implementation in the 
computational procedure are given in Sec. 3; Sec. 4 illus- 
trates the 4 types of chemical structures that were included, 
by displaying one of each type, all the others being collected 
in Appendix B; Sec. 5 assembles the corresponding spectra 
in 4 graphs, one for each type, and compares the resulting 
overall emission spectrum with an observed galactic spec- 
trum. Section 6 discusses the defects of this simulation and 
Sec. 7 explores possible improvements. 



2 RADIATED ENERGY 

Present chemistry simulation codes do not include sponta- 
neous radiation. A simple analytical argument can, however, 
be made to deliver an operational approximation to reality. 
When energy is deposited in a molecule, several modes are 
excited; which ones depends on the excitation process. If the 
modes were strictly harmonic and, therefore, completely un- 
coupled, each and everyone would start radiating in the IR, 
at its own rate, different than the others. Assuming the exci- 
tation rate to be slow enough, each mode would ultimately 
radiate the energy which was initially deposited in it; the 
radiated power would be obtained by multiplying this with 
the frequency of excitation events. But, in a real molecule 
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above its ground state, anharmonicity cannot be overlooked 
and the modes are coupled more or less. As a consequence, 
intramolecular vibration redistribution then quickly changes 
the initial distribution into a new one, which is mostly inde- 
pendent of the excitation process, but characteristic of the 
molecule itself when in dynamical equilibrium (see Papoular 
2001, 2002). 

Internal relaxation strongly impacts the spectral emis- 
sion, as we now show on a simplified ideal case. Assume first, 
for simplicity, that the molecule has only 2 vibration modes, 
and let Ei be the energy, n the radiative lifetime of the ith 
mode, with i = 1,2, and e the energy exchanged between 
modes. Also let c" 1 represent the fraction of E e xc (the total 
excitation energy) that is stored in mode i when internal 
relaxation is completed, i.e. the above mentioned character- 
istic energy distribution in dynamic equilibrium. Then, 



Ei = -Ei/n+e 
E 2 — —E 2 /r 2 — e 
e = a(c 2 E 2 — ciEi), 



(i) 

(2) 
(3) 



where the dot designates the derivative and a is the 
mode coupling constant. Note that, when Ei oc c^ 1 , energy 
exchange stops. Internal and radiative relaxation may be 
represented by a solution to this set, of the form 



Ei = aiexp(—t/r) + f3iexp(— t/cr) 
E 2 = a 2 exp(—t/r) + f3 2 exp(—t/a) 
e = aexp(—t/r) + f3exp(—t/a). 



(4) 
(5) 
(6) 



The meaning of these equations is clear when, as is usu- 
ally the case, the coupling constant, a, is strong enough that 
a <C r. In a first, short, intramolecular relaxation phase (the 
transient, with time constant a), radiation is insignificant, 
but the energies in the two modes are forced to quickly read- 
just (from Eoi to ~ on) so as to nearly annul e (eq. 3). They 
then can relax at the same, much slower, rate, r, intermedi- 
ate between the two radiation rates n and r 2 , while energy 
is being continuously transfered from the less radiative to 
the more radiative mode. 

That a < r, results from the fact that the radiative 
lifetime in the IR is very long (>lms), while the intramolec- 
ular vibrational relaxation is shorter than 1 ns, in general, 
as illustrated by numerical modeling on several molecules 
(see Papoular 2001, 2002 and 2006). 

The 8 unknown constants appearing in the solution eq. 
4 to 6 are determined by first substituting in eq. 1 to 3, and 
setting to the sum of the coefficients of the exponentials 
in a and a, separately and for each equation. Then, initial 
conditions are applied: Ei(t = 0) = Eqi. The results are 
cumbersome to an unwarranted degree unless d C r. In 
that case, the procedure leads to 



(7) 



Moreover, the intra- molecular and radiative phases 
hardly overlap, so that, at the start of the latter, there is 
no more energy exchange between modes: e = 0, so that 

c\ol\ — c 2 a 2 . (8) 

Then, the required constants are approximately given 

by 



oi — c 2 - 



Eoi + Eo 2 

ci + c 2 
Eqi + Eq 2 



-132 = 13: 



(9) 

Q2 = C] ; (10) 

(11) 

(12) 
(13) 

a = a±a 2 . (14) 

OL 2 T\ + Ct\T 2 

All of the above is borne out by numerical resolution of 
system 1 to 3. Finally, the energies radiated by modes 1 and 
2 (after extinction of the transient) are, respectively 



ci + c 2 

C\Eq\ — C 2 Eq 2 
Cl + C 2 
_ T\T 2 (C\ + C2_) 

cm + c 2 r 2 
1 



o(ci + c 2 ) 

Tl - T 2 



J Eidt/ri = roti/ri 
J E 2 dt/r 2 — ra 2 /r 2 = 



ci + c 2 
cm + c 2 r 2 

Cl + C2 



-T 2 OLl 



cm + c 2 r 2 



ria 2 , 



(15) 
(16) 



(17) 



whose sum equals ai + a 2 , as it should, since there is no 
outlet for the initial energies other than radiation. The ratio 
of total energies radiated at the respective mode frequencies 
is 

Oil t 2 _ c 2 r 2 
o 2 n ci n 

The forms of the solutions above are easily generalized 
to any number of modes, assuming a uniform coupling con- 
stant. This is inocuous since the coupling constant only im- 
pacts the transient time constant, cr, which remains very 
short anyway (a few tens of picoseconds). 

The radiative lifetime is the inverse of Einstein's spon- 
taneous emission coefficient (see Wilson et al. 1955) 

-m(is)\ (18) 



w 64vrV , x2 



3/ic 3 

in electrostatic CGS units, and written in simplified form, 
for one particular transition and overlooking degeneracies. 
Here, m(y) is the dipole moment oscillating component at 
frequency v. Now, it comes as no surprise (see Atkins 1997) 
that the corresponding IR intensity is also linked to m(y) 
by 

I[y) = 2&vm(yf , (19) 

where m is in Debye and v in cm -1 . Thus, Einstein's co- 
efficient scales like v 2 l(y) and is 3 m(is) 2 . Finally, the energy 
ultimately radiated by mode v scales like 



E* oc is 3 c 1 (is)Eex C m(is) 2 



(20) 



where the mode index i was replaced by its frequency. 
Equation 20 can also be derived from the classical oscilla- 
tor emission by multiplying with the oscillator strength (see 
Kuhn 1962, Wilson et al. 1965). Its range of application is 
therefore wider than those based on Einstein's coefficient or 
integrated line absorption, as these are restricted to transi- 
tions from or to the ground state, and they are delta func- 
tions unaffected by anharmonicity, so they cannot account 
for line broadening and shift. 
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Appendix A below illustrates this procedure as applied 
to ethylene C2H4, whose sparse IR spectrum makes for a 
better presentation of the effects of anharmonicity, and dif- 
ferences between IR emission, absorption and energy distri- 
bution. 



3 THE COMPUTATIONAL PROCEDURE 

The present work is based on the use of various algorithms of 
computational organic chemistry, as embodied in the Hyper- 
chem software provided by Hypercube, Inc., and described 
in detail in their publication HC50-00-03-00, and cursorily, 
for astrophysical purposes, in Papoular (2001). Here I use 
the improved version Hyper 7.5. While almost all the com- 
putations used standard algorithms, some required writing 
special purpose subroutines, a capability also provided by 
the software. Semi-empirical methods (PM3 and AMI) were 
used for molecules up to about 100 atoms in size; their accu- 
racy for spectral frequencies is of a few percent. Beyond that, 
the less accurate, but faster Molecular Mechanics (MM+) 
method was preferred. Anharmonicity is built in and tailored 
by comparison with laboratory experiments on hydrocarbon 
molecules. 

Once a molecule has been selected, the first step is the 
building and optimization of its structure, i.e. minimizing 
its potential energy by automatic standard procedure, at 
K (no kinetic energy in the molecule). Next, a type of exci- 
tation is adopted. An isolated molecule in space can acquire 
energy by, for instance, absorbing a photon, or by one of its 
dangling (C-) bonds capturing an H atom, thus increasing 
its potential energy by an amount equal to the binding en- 
ergy (Papoular 2012). In the following, we assume that the 
average interval between two successive excitation events is 
much longer than the IR lifetime of the molecule. So, the ra- 
diative power of a molecule is the product of the frequency 
of exciting events and the energy radiated by the molecule 
upon radiative relaxation. Although different excitation pro- 
cesses were studied, this paper focuses on H capture by one 
of the dangling C bonds that are likely to be present at 
any time at the molecular periphery. In this process, a free 
H atom approaches the grain closely enough to feel the at- 
traction of a carbon atom whose bonds are incompletely 
occupied. Under favorable initial conditions, it falls into the 
corresponding potential well. It was shown that this cap- 
ture reorganizes the distribution of the constitutive atoms 
in space so as to lower the ground state of the newly com- 
pleted structure by an amount equal to the CH bond energy 
(Eb ^4.5 eV). This is equivalent to depositing this energy 
in this new structure, an energy which is immediately and 
evenly divided into kinetic energy and potential energy (rel- 
ative to the ground state of the completed structure). The 
kinetic energy itself, initially localized in the CH stretching 
vibration, is then redistributed among this and other (not 
necessarily all) molecular vibrations. This could not occur 
were it not for the couplings between modes. While normal 
modes are uncoupled by definition, the vibration modes be- 
come coupled as soon as some energy is deposited in the 
molecule, because of anharmonicity, which sets in above the 
ground state. Note that, in this process, by contrast with 
UV photon excitation, the molecule remains in the electronic 
ground state even though atomic vibrations are excited. 



To simulate this process on the computer, we start with 
an optimized structure at rest, with all bonds occupied. One 
of the CH bonds is then rotated out of the molecular plane 
by several degrees and extended from its original length by a 
few hundredths of A so as to increase the molecular poten- 
tial energy by, say, E exc , which is tantamount to depositing 
this energy in the molecule. The molecule is then set free to 
perform a molecular dynamics run in vacuum, during which 
all atomic motions are monitored on the computer screen 
and recorded, as well as the total kinetic and potential ener- 
gies, as a function of time. This first run thus simulates the 
first, transient, phase described by eq. 4 to 6. 

Initially, as expected, the dominant motion is the CH 
stretch accompanied by excursions out of the molecular 
plane. Soon, however, the other atoms are dragged along: 
inspection and comparison of the variations of, say, the bond 
lengths immediately reveals that energy constantly flows 
from bond to bond, thus increasing the amount of energy 
shared by other bonds. This process translates into a con- 
tinuously decreasing standard deviation, a, of each and ev- 
ery variable, like the total kinetic energy of the system, for 
instance. After some time (the relaxation time), a dynamic 
equilibrium sets in, which is characterized by constant aver- 
age relative standard deviations. The relative standard de- 
viation of the energy, cfe/E is roughly N^ 1 ^ 2 , where Nm 
is the number of effectively excited normal modes (see Reif 
1984). The relaxation time is highly variable, depending on 
the grain structure and size, and on its initial perturbation. 
It may be as long as hundreds of picoseconds. This is still 
much shorter than the minimum radiative lifetime (~ 1 ms). 
The energy radiated during this time is negligible. 

At the start of this relaxation phase, the initial energy 
distribution among modes depends mainly on the excitation 
mechanism, but it tends asymptotically toward a different 
final distribution (the c's). I call this the "characteristic" 
distribution for it depends mainly on the molecular structure 
and composition. 

In dynamic equilibrium, energy is still continuously ex- 
changed between modes, but is essentially conserved in the 
molecule as radiation is not included in the PM3 code. The 
time averaged energy distribution can be identified with the 
c _1 's of Sec. 2. Once this state has been reached, a new run 
is launched, to simulate the emission phase, during which 
all variables of interest are monitored. First, we focus on 
the instantaneous total kinetic energy. The FT of its excur- 
sions from its mean delivers the excitation energy spectrum 
c _1 (y)E eX c required in eq. 20. It exhibits peaks in the vicin- 
ity of the normal modes. Over and above the physical broad- 
ening due to anharmonicity, these peaks are also artificially 
broadened by the finite length, A t, of the run: this contribu- 
tion to the overall width at half maximum (FWHM, cm" 1 ) 
is rsj (cAt) -1 . The longer the run, the less significant will this 
contribution be compared with the physical broadening: in 
general, 10 ps will do. 

The dipole moment of the molecule is similarly mon- 
itored in dynamic equilibrium, and its Fourier transform 
similarly gives the dipole moment spectrum , ra(z/), also re- 
quired in eq. 20, which can finally be computed to give the 
desired emission spectrum. 
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Figure 1. One member of the "aromatics" class of selected struc- 
tures. Atoms are not labeled, but C and H atoms are indicated by 
black and green bond segments, respectively. Aliphatic chains are 
expected to naturally bond to pure aromatics, and help slightly 
shift the spetral lines of the latter. Other members in Fig. 18. 




Figure 4. One member of the "trio" class. Other members in Fig 
22 to 24. 
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Figure 2. One member of the "aliphatics" class. Another member 
in Fig. 19. 



4 THE SELECTED MODEL UIB CARRIERS 

The attempt at simulating UIB spectra, as described in the 
present work, required the computation of emission from 
a large number of candidate structures so as to determine 
the characters of their possible contribution to UIB emis- 
sion (see Papoular 2010). From this data base, 21 different 
molecules were selected, falling in 4 categories: aromatics, 
aliphatic chains, oxygen-bridged chains and "trios" (essen- 
tially a 5-membered ring squeezed between two 6-membered 
rings) . These denominations are used here in a very restric- 
tive sense and apply only to the present, or slightly modified, 
structures. In this section, one member of each carrier cat- 
egory is illustrated (Fig. 1 to 4), The remaining structures 
are illustrated in a Appendix B (Fig. 18 to 24). 

The structure illustrated in Fig. 21 differs fundamen- 
tally from the others in that it is a cluster of four O-bridged 




Figure 5. All 5 aromatic emission spectra. 

chains (Fig. 3 and 20) bonded together by van der Waals 
forces. Its IR spectrum does not differ essentially from the 
sum of the spectra of its components. This suggests that the 
model UIB carriers proposed here could be taken separately 
or in clusters. 



5 THE EMISSION SPECTRA 

Figure 5 to 8 each show the emission spectra of all members 
of one class, using one color for each member. The weighted 
spectra of the 21 selected molecules were added to give the 
expected overall emission. The weights were then tailored 
to steer the sum towards a typical observed and published 
galactic spectrum chosen for its quality and spectral extent 
(galaxy NGC1482; Smith et al. 2007). Figure 9 displays the 
result, together with the observed spectrum for comparison. 
Figure 10 shows the effect of smoothing the spectrum of 
Fig. 9. This simulates the expected effect of increasing the 
number of similar, but slightly different structures, of the 
selected set. 



Figure 3. One member of the "O-bridged chains". Oxygen is 
colored in red. Other members in Fig. 20 and 21. 



6 DISCUSSION 



It must be stressed again that the emission spectra displayed 
here are naturally broadened by the inherent anharmonicity 
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Table 1. Atomic composition of dust 
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of the molecular dynamics and the ensuing coupling between 
atoms, as described by the computational codes used in this 
work. The only artificial broadening involved in the calcu- 
lation is due to the finite length in time of the molecular 
dynamics run, from which the dipole moment spectrum is 
deduced by FFT. For a 10 ps run, the broadening is 0.0035 
/mi, which does not impact the spectra notably. 

When assessing the UIB simulation, it must be remem- 
bered that the inevitable defects of the chemical simulation 
codes entail an error of a few percent on the computed fre- 
quencies. With this in mind, it appears that the synthetic 
model spectrum exhibits a limited number of bands falling 
respectively at 3.25, 5.5, 5.65, 6.3, 7.4, 9, 11.4,12.8, 17.1 and 
18.1 /mi, near the main generic UIBs: (3.3,) 5.3, 5.7, 6.2, 
7.7, 8.6, 11.3, 12.7 and 17-18 /mi. The largest discrepancies 
occur for the 7.7 and 8.6 /mi bands, and are of order 5 %; 
they overemphasize the dip in between. 

In Fig. 11, each colored line represents the sum of the 
weighted spectral contributions of the members of one class 
of structures. With the exception of the CH-stretching peak, 
each peak appears with a different color. This confirms that 



each of the 4 types of structures mainly contributes to a 
different UIB band, or set of bands: 

-aromatics: 11.3, 17-18 /mi 
-aliphatics: 8.6 /mi 
-O-bridged chains: 7.7 /mi 
-trios: 5.3-5.7, 6.2, 12.7 /mi. 

These correlations stem from specific atomic motions 
favored by the structures. More detailed assignments are 
given in Papoular (|2010p and were used to select the 21 
structures considered here. These findings help tailoring new 
models to fit other observed spectra, which may differ in 
their relative band intensities. An interesting case in point 
is M82, in the halo of which, the ratios of the 3.3, 3.4 and 6-9 
m um bands were found to vary significantly with distance 
from the Center (see Yamagishi et al. 2012). 

Table 1 gives the total number of atoms, the number of 
C, H, O and S atoms in each molecule, and its weight, w, in 
the sum. The last row shows that, for this combination of 
molecules, 

[O]/[C]=0.13, [H]/[C]=2, [S]/[H]=0.0015. 
Not unexpectedly, these numbers correspond, in Van 
Krevelen's diagram, to kerogens type I (see Durand 1980, 
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Figure 6. The 2 aliphatic emission spectra. 
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Figure 9. Black: the overall emission spectrum, obtained by 
adding all 21 member spectra, weighted by the coefficient w, tab- 
ulated in the last column of Table 1; area under curve is filled. 
Red curve: the flux curve, Lu , for NGC1482, from Smith et al. 
(2007); magenta dots mark the feature peaks tabulated in their 
Tab. 3; the continuum is not subtracted. 
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Figure 7. All 6 O-bridged-chain emission spectra. 
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Figure 8. All 8 trio emission spectra. 



Figure 10. Black: FFT smoothing of the same overall emission 
spectrum as in Fig. 9; resolution: 0.1 /im . Red curve: the flux 
curve, I u , for NGC1482, from Smith et al. (|2007|> : magenta dots 
mark the feature peaks tabulated in their Tab. 3; the continuum 
is not subtracted. 



pi 22- 126). Such kerogens are the least deeply buried in earth 
and, hence, the least aromatized (or mature, or "evolved" 
or "coalified") ones. Structural models for these, and other 
kerogens, were developed long ago; as stated in the Intro- 
duction, these models helped select the emitter structures 
considered in the present work. 

Structural models developed for kerogens and coals in- 
clude of course many more types of structures. For our 
present purposes, the choice is further constrained by as- 
tronomical data. Thus, alkene chains were not included, as 
double CC bonds carry strong features at 5.27 and 5.7 /xm, 
whereas the corresponding features observed in space are 




Figure 11. Breaking the overall emission spectrum by classes. 
Red: aromatics; green: aliphatics; blue: O-bridged chains; black: 
trios. 



Figure 12. The 3.3 /im region. 



relatively quite weak. Similarly, cyclohexanes (saturated 6- 
membered rings) could not be put to use to any extent. 

Note that nitrogen and OH tails were used abundantly 
in modeling the 21- and 30-35-/xm features of protoplanetary 
nebulae (Papoular, 2011). They were not used here, as these 
features do not show up in the spectrum of NGC1482 pub- 
lished by Smith et al. (2007). These authors also indicate 
that the conspicuous 18.9-/xm band is mainly due to the ion 
SIII, which justifies its absence from the synthetic spectrum 
of Fig. 9. 

One of the most conspicuous discrepancies between the 
synthetic and observed spectra is that the envelopes of the 
synthetic (black) "bands" in Fig. 9 are not as smooth as 
those of NGC1482 (red). That is made clearer by zooming 
on the main bands, as in Fig. 12 to 15. Obviously, this is 
due to the sparseness of the selected model structures: the 
total number of atoms in the 21 selected molecules is 1068, 
giving rise to only 3198 vibration modes. It is clear that very 
many more structures can be built along the lines suggested 
by Fig. 18 to 24, to try and match the infinite diversity to 
be expected in space. But this must be left for later work. 
A glimpse of the effect this would have on the spectrum of 
Fig. 9 is gained by smoothing, as in Fig. 10. 

This is also a reminder of the fact that other selections 
of model structures may well fit observed spectra: in the 
present model, they are probably all present in space. 

The close up views of the model bands reveal that the 
bunching of many peaks within a given band is necessar- 
ily accompanied by the appearance of underlying plateaus, 
formed by superposition of their wings. Also note that, be- 
tween 6 and 15 /im, the sum spectrum never falls to zero, 
meaning that the selected set of molecules contribute to the 
underlying continuum. According to this model, therefore, 
the underlying continuum due to other families of carriers 
must not necessarily be drawn through the troughs between 
bands. 

Another defect of the model is the excessive relative 
intensity of the 3.3-/xm band: it is about 5 times more in- 
tense than the 7. 7- /mi band. Although it is not included in 




Figure 13. The 5.5-and 6.2-/mi region. 



the NGC1482 spectrum for comparison, other observations 
suggest that the intensity of this band is of the same order 
as the 6.2- and 7. 7- /mi bands (see Peeters et al. 2004b). In 
the model spectrum, the main contribution to this band is 
that of the trio class, suggesting that this could be reme- 
died by replacing some of the H terminations by methyls 
(or CH 2 chains). A new structure was therefore built from 
trio a (Fig. 4), by replacing all the H's by CHs's (except 
those which cap the two pentagons) . Their emission spectra 
are compared in Fig. 15: as expected, the CH stretchings are 
replaced by symmetric (umbrella) stretching of the methyls 
(to the left, probably because of code inaccuracy) and , to 
the right, symmetric stretchings of the two CH's capping a 
pentagon, and asymmetric stretchings of methyls. For co- 
existing H and CH3 terminations, these aliphatic features 
could conveniently account for the 3.4-/i massif as well as 
the underlying plateau observed in the sky. Further mod- 
eling requires more simultaneous observational data on the 
near- and mid-IR domains. 
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Figure 14. The 7.7-and 8.6-/mi region. 
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Figure 15. The 11.3- , 12.7- and 18- /im region. 

With one exception (Fig. 19), the carriers considered 
here are all chemically bonded, medium sized structures. 
Molecules in space must be expected to suffer mild encoun- 
ters and coalesce under Van der Waals attractive forces, to 
give clusters like in Fig. 19, which may become progenitors 
of dust grains. The study of such clusters should help assess 
the relative contributions to UIBs, of clusters (grains) and 
free molecules. 



7 CONCLUSION 

Clearly, more work is needed to obtain a satisfactory fit of 
model to observed spectra. The straightforward part might 
consist in including many more structures similar to those of 
the 4 mentioned classes. Exploring larger structures would 
help, as diversity increases with size, but this is more chal- 
lenging. Also, that would call for more accurate and efficient 
codes and machines. 

Studying the transition from molecules to grains, how- 



Figure 16. The effect of replacing H terminations of trio a 
(mauve) by methyl terminations (green): the "3.3"-/im band is 
replaced by much weaker massifs on its wings, and the "6.2"-/im 
band is considerably increased. 

ever, cannot be avoided if important issues are to be ad- 
dressed, such as the emission continuum, the graphitization 
at the origin of the 2175 A extinction band and the evolu- 
tion from the aromatic 3.3-/im to the aliphatic 3.4-/xm band. 

Even in the realm of medium-sized molecules, it should 
be possible to consider more tightly knitted structures than 
those studied here, which are essentially 2D. By contrast, 
taking on ionized molecules at present seems very laborious 
and uncertain with available simulations codes for molecular 
dynamics. 

From the strictly chemical point of view, it would be in- 
teresting to understand why and under what circumstances, 
structures similar to those described in this work would come 
to be privileged. 

Laboratory experiment could help advance these issues. 
In particular, the IR emission of a gaseous assembly of 
medium-sized molecules could be collected and spectrally 
analyzed at different pressures and temperatures: this may 
provide both a test of emission theories and an understand- 
ing of molecule-to- grain transition. 



8 APPENDIX A: LINE BROADENING 

Figure 17 is an example of spectral results obtained by 
the procedure outlined in Sec. 2. It is a zoom on the CH 
stretching band of the ethylene molecule after excitation by 
bond extension (which injects 2.5 Kcal/mol in the molecule). 
The bottom graph represents the spectrum of the energy or 
power radiated by the molecule. This was obtained by mon- 
itoring the total electric dipole moment of the molecule and 
applying eq. 20. For comparison, the middle graph repre- 
sents the normal modes of ethylene, and the upper curve 
is the FT of the excursions of the kinetic energy from its 
mean, and represents the spectral distribution of energy in 
the excited molecule. Most normal mode lines are conspic- 
uously absent in emission , because the corresponding vi- 
bration modes are either IR-inactive or not excited. As a 
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Figure 19. The other member of the aliphatics family. 



Figure 17. The CH stretching band of ethylene, obtained by 
the procedure outlined in Sec. 2. From bottom to top: radiated 
energy (or emission); normal modes (only 2 are IR active); energy 
distribution among modes. 




Figure 18. Other members of the aromatics family. 

consequence of excitation and anharmonicity, the separa- 
tion between the two IR-active CH-stretching lines is nearly 
twice that of the corresponding normal modes, and the 
shorter wavelength emission peak is dominant. The hori- 
zontal segment, of length 0.0012 /im (bottom graph), repre- 
sents the artificial broadening due to the finite calculation 
run time (30 ps). The real broadening, due to anharmonic- 
ity and mode coupling, is ^0.004 /mi. The relative shift of 
~ 0.035/im (1 %) between the middle spectrum and the 
other two is probably to be ascribed to minor inconsistencies 
between the two distinct algorithms which compute normal 
modes and molecular dynamics, respectively. 




Figure 20. Other members of the O-bridged chains family. 
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9 APPENDIX B: SELECTED MOLECULES 

Figures 18 to 24 display all the chemical structures used used 
in the present work to model the spectrum of NGC 1482, 
except the head members of each class, displayed in Fig. 1 
to 4. 
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Figure 21. This last member of the O-bridged chains family is a 
cluster of chemically bonded molecules, held together by van der 
Waals forces. 
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Figure 22. Other members of the trios family. 




Figure 23. Other members of the trios family. 
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Figure 24. Other members of the trios family. 



